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Abstract 

Continuum Monte-Carlo simulations at constant pressure are performed on 
short chain molecules at surfaces. The rodlike chains, consisting of seven 
effective monomers, are attached at one end to a flat twodimensional sub- 
strate. It is found that the model exhibits phases similar to the liquid con- 
densed and liquid expanded phases of Langmuir monolayers. The model is 
investigated here for a wide range of pressures and temperatures using a 
special form of constant pressure simulation compatible with the symmetry 
breaking during tilting transitions in the liquid condensed phases. At low 
pressures the chains undergo a tilting transition exhibiting tilt directions 
towards nearest and also next nearest neighbours depending on tempera- 
ture. At elevated temperatures and low pressure the film enters a fluidlike 
phase similar to the liquid expanded phase observed in experiment. 

PACS: 68.10.-m, 68.16.+e 



I. INTRODUCTION 



Amphiphilic molecules, such as fatty acids {CH^{CH2)nC00H) or phospholipids, 
form monomolecular layers on water surfaces in which the hydrophobic headgroups are 
immersed into water while the hydrophobic alkane chains point outwards [|I|. This phe- 
nomenon is utilized in the Langmuir-Blodgett coating technique, in which under suitable 
conditions on pH and temperature, monolayer or multilayer films can be transferred to a 
solid substrate by successively dipping the substrate into a compressed layer on a water 
surface. Langmuir-Blodgett films are highly ordered and their molecular ordering depends 
on the structure of the compressed layer on the water surface prior to dipping f^. The 
ordering of Langmuir-Blodgett films plays an important role in applications ranging from 
microlithography and electrooptics to biochemical sensing and tribology. 

Many experiments have in recent years investigated the molecular structure and phase 
behaviour of Langmuir monolayers (see [Q-]^ for reviews). These studies were spurred 
by the possibility of combining modern x-ray diffraction methods with more traditional 
thermodynamic measurements 0,0] • The rodlike non axially symmetric molecules give 
rise to a large variety of phases and phase transitions |Q with a rich polymorphism. 

Despite differences in detail monolayers of fatty acids, phopholipids, alcohols and 
esters exhibit a similar phase diagram This fact alone suggests that the structural 

similarities between amphiphilic molecules are responsible for the similarities observed 
in the phase diagrams. Hence the chemical and atomistic differences between different 
molecules may not be necessary to understand the phase diagram of amphiphiles. This 
observation motivates our studies of highly idealized coarse grained models to elucidate the 
mechanisms underlying various transitions. Among these transitions the tilting transition 
in the liquid condensed phases has been our main interest. We study idealized models 
in order to separate the influence of various coarse parameters (such as chain flexibility. 
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symmetry of the molecules, mobility and geometry of head groups etc.) from that of finer 
details (such as interaction potentials, cis-trans conformations, hydrogen atoms etc.) on 
the tilting transition in these monolayers. 

Given the advances in numerical simulation techniques there have already been numer- 
ous simulation studies of chemically realistic P^TP[ and idealized coarse grained models 
p0|^5| for the tilting transition in amphiphilic monolayers. Chemically realistic models 
with much atomistic detail have an advantage in their ability to describe dense packing 
effects more realistically. Molecular dynamics simulations of such models are very useful 
for elucidating the behaviour of such models away from phase transitions. Near a phase 
transition Monte Carlo methods can be more efficient for obtaining the statistical aver- 
ages. Highly idealized coarse grained models may be simpler (and sometimes faster) to 
simulate but they often suffer from neglecting some aspects of the physics. For example, 
treating the alkane chains as rigid cylindrical rods, originally suggested in ||2^ , disregards 



the melting transition of the layer, and the intrachain conformational disorder. Such rigid 
rods grafted onto a triangular lattice were studied extensively in pT|J23J2^ using contin- 



uum Monte Carlo calculations. While this model exhibits a tilting transition it cannot 
show a liquid expanded phase or the restructuring and melting of the head group lattice. 
Although the grafting of the head groups in this model severely limits the applicability 
of the model to Langmuir layers the model is useful and realistic for chains that are per- 
manently grafted to a solid substrate [0. Chemically realistic as well as coarse grained 
models suffer from the general limitation that the effective potentials are only inaccurately 



known 26 . 



Layers of flexible chain molecules with mobile head groups p7| , p8| are intermediate 
between chemically realistic united atom models and idealized coarse grained models 
such as rigid rods or even lattice models Such models were studied in using 



continuum Monte Carlo simulations at fixed volume. It was found that the condition of 
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fixed volume imposes a symmetry breaking field which can suppress the restructuring of 
the head group lattice during the tilting transition. This was concluded from the presence 



of metastable boundary induced defect structures While such defects are known to 
exist in experiment and are therefore of interest in their own right [^], they make the 
extraction of thermodynamic equilibrium properties from small scale simulations very 
difficult. Hence it is of great interest to perform also simulations at constant pressures in 
which the symmetry breaking field is absent. 

In the present paper we report simulations at constant pressure using an algorithm 
which is compatible with the symmetry breaking during tilting transitions. To facilitate 
comparison we investigate a model similar to the one studied at constant area in p7|j28| 



Our main objective is twofold: Firstly we would like to confirm the existence of an ex- 
pected orthorhombic distortion of the head group lattice during a tilting transition in a 
coarse grained model. Such observations are essential for the interpretation of constant 
volume simulations in which the distortion is suppressed. For example, the impossibil- 
ity to distort or contract the simulation box in such simulations combined with periodic 
boundary conditions can lead to the appearance of negative pressures at low densities if 
the interaction potential has an attractive part. Secondly, we want to understand the se- 
lection mechanisms for nearest versus next nearest neighbour tilt direction and the rather 
high transition temperature of the first order melting transition observed in previous con- 



stant volume simulations of the same model p7| , p8| , |31| , p9| . In [^3)111 fiuidized phase 
appeared at temperatures roughly four times higher than the tilting transition tempera- 
ture. This seems to rule out the possibility of identifying the molten state with the liquid 
expanded phase of Langmuir monolayers. An important result of the present study is to 
show that in constant pressure simulations the fiuidized phase appears already at much 
lower temperatures. 
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II. MODEL DESCRIPTION 



Each amphiphilic molecule is represented through a chain of seven effective monomers 
labeled i — 0, ...,6 with i — corresponding to the hydrophilic head group of the am- 
phiphile. The cartesian coordinate system in threedimensional space is chosen such that 
the headgroup i = is restricted to move in the z — plane representing the twodimen- 
sional substrate. All the effective monomers are connected through a finitely extendable 
nonlinear elastic (FENE) potential 



where d is the distance between monomers, and > is the spring constant. The FENE 
potential is harmonic at its minimum do but the bonds cannot be stretched beyond a 
maximum length determined by di,i. The stiffness of the rodlike molecules whose alkane 
chains are predominantly in an all trans conformation is incorporated by a bond angle 
potential 



where Cba is the force constant and 6i is the angle formed by the three monomers 
Note that Vba is a three body interaction. All monomers except nearest neighbours within 
the same chain interact through a Lennard- Jones potential. The Lennard- Jones potential 
is truncated and shifted such that it vanishes at the truncation point. If e is the interaction 
strength and a its range then 






(2) 



e{{a/dY' - 2{a/df - {l/dLjf^ + 2{l/ dLjf) 



for d < dijcr 



VLj{d) = 



(3) 







for d > dijcr 
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where d = dijcr is the truncation point. 

The model above is similar to the model used for the constant volume simulations in 
27| , |28[] . It extends and generalizes an earlier study of a system consisting of perfectly rigid 



rods grafted to a hexagonal lattice and interacting with Lennard- Jones interactions ||2 
The main difference with [ P7|P5 | consists in simulating at constant spreading pressure 



rather than constant substrate area, and in replacing the cutoff harmonic potential for 
the bondlengths with a FENE potential. 

The earlier simulations at constant volume were found to lead to boundary induced 



defects ||29|. These can significantly perturb the long range correlations at phase transi- 
tions. Most molecular dynamics simulations are similarly carried out at constant volume 
or aspect ratio of the simulation box. Simulations at constant pressure on the other hand 
reproduce the experimental situation more exactly, and allow the defects to relax via 
distortions of the head group lattice. 



III. RELATION BETWEEN THE MODEL AND EXPERIMENT 

The model presented in the previous section is not intended to be a chemically realistic 
model for Langmuir monolayers. On the contrary we wish to study a highly idealized 
coarse grained model in order to understand the influence of varying degrees of idealization 
on the mechanisms responsible for the tilting transitions in the condensed phases. The 
model has resulted from previous systematic investigations of even more idealized models 
21|,|23|-|5],|^|8|. 



Although the model is not intended to be chemically realistic it is of interest to connect 
the model with reality through a rough mapping of length and energy scales. Such a rough 
mapping can be helpful in deciding whether the phase transitions of the model may be 
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related to the phase transitions of the experimental system. 

If a typical fatty acid with chain lengths from 12 to 30 carbons is represented by a 
model chain with seven effective monomers then each effective monomer of our model 
chains must represent roughly between two and five methyl groups. Hence, in such a 
mapping the diameter of the effective monomers is of order lOA, and the Lennard- Jones 
interaction range a is of order 5A. For the energy scale a value for e of order 200K is a 
very rough estimate. 

The simple attachment of the head groups to the substrate represents an important 
overidealization of our model. Experimentally the form and interactions of the head 
groups appear to play an important role. While our model is more realistic than fully 
discrete lattice models ||2^ it contains less chemical detail than united atom models which 



have been investigated using molecular dynamics simulations |^,|T^,0. The main differ 



ence is that our molecules are axially symmetric while the zigzag structure of united atom 
models destroys this symmetry. Consequently we do not expect the model to reproduce 
phases associated with a herringbone ordering. We feel, however, that our model while 
neglecting these features represents a good compromise between computational efficiency 
and chemical realism. 



IV. SIMULATION DETAILS 



Lengths and energies in our simulations will be measured in dimensionless units defined 
by setting cr = 1 and e = 1. In these units the other parameters for the simulations are 
chosen as do = 0.7, du = 0.2, dij = 2.0 and cu = 100, Cba = 10 which ensures that the 
chains do not intersect. The choice Cbi,Cba ^ e produces stiff rodlike chains. The system, 
consisting of 144 chain molecules, is confined to an area of sidelengths L^, Ly in the z = 
plane. The height of the simulation box ^ Gd^ is chosen much larger than the length 
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of the chain molecules. 



A simplified groundstate analysis of the model was performed previously p7|j28|j3T 



It predicts a tilting transition between a tilted and an untilted state, and highlights the 
importance of packing effects. 

The continuum model is simulated in a canonical ensemble in which temperature, 
spreading pressure and particle number are kept constant. The simulation is carried out 
using a Metropolis Monte Carlo procedure in which individual monomer positions and the 
area A = L^Ly of the simulation box are updated in continuous space. The continuous 
position space requires to use methods adapted from molecular dynamics simulations 
for evaluating the interaction energies. The maximal jump distance of monomers in a 
single move is chosen to optimize the acceptance rate while ensuring that the chains 
cannot intersect each other. The cutoff in the interaction potentials allows us to use an 
adaptation of the link cell algorithm to keep track of interacting neighbouring monomers. 
The simulation box is subdivided into cells of sidelength larger than or equal to the 
interaction radius. Particles in a subcell can only interact with particles in neighbouring 
cells. The cell contents are stored using specially designed linked pointer lists. The 
algorithm was developed and tested for constant volume simulations [^,^,^, and its 
speed depends strongly on the density and the interaction ranges. 

The partition function Z^vn^T of the isobaric ensemble is obtained from the isochoric 
partition function Z^at through Laplace transformation. Denoting the surface area by 
A and the spreading pressure by 11^ it reads 

/•oo r 

Znu^t = J^ J exp{-f3{UAA + V{r))) drdA (4) 

where (3 = l/{kBT), is the Boltzmann constant, dr is the integral of all coordinates, 
and all constant prefactors have been suppressed. Expectation values for an observable 
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X{r) in this ensemble are calculated as 

{X) = — r [ X(s) exp (-/3(V(s) + UaA) + N\nA) ds dA (5) 

in terms of the scaled coordinates s = {tx/Lx, Ty/Ly, Vz) of all the particles. The expecta- 
tion values can be calculated with the usual Metropolis scheme by introducing an effective 
Hamiltonian 

i/e//(s, Ly) = V{s) + TlALxLy - NkBTln{LxLy) (6) 

containing two additional degrees of freedom and Ly. Attempts to change the area 
A = L^Ly of the simulation box are carried out after each Monte Carlo step (IMCS = 1 
update/monomer) of regular coordinate updates. A maximal step size of AL^ — A.Ly — 
±0.01 for the area moves was found to yield optimal acceptance rates for the chosen 
system size. The evaluation of the energy difference AV{s) for an area move involves 
a recalculation of all particle interaction potentials because the coordinates are rescaled 
differently in the x- and y-directions. This makes area moves computationally expensive 
and leads to a deceleration of the constant pressure algorithm by a factor of roughly two 
compared with constant area simulations. 

The simulations are started from an untilted configuration of 12 x 12 chain molecules. 
The chains are positioned on a triangular lattice and directed perpendicular to the sub- 
strate. If a tilted initial configuration with hexagonal headgroup lattice is used at low 
temperatures the system was found to evolve into metastable states with lifetimes well 
beyond 10^ Monte Carlo steps. In all simulations periodic boundary conditions are ap- 
plied in the xy-directions. In each run the system was first equilibrated for 20000 Monte 
Carlo steps (updates per monomer). Subsequently averages were recorded every 500th 
MCS over a period of 50000 MCS. These calculations consumed several 100 hours of CPU 
time on IBM RS6000 370 equipment. 
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The two dimensional pressure tensor was estimated form the forces on each particle 
using the virial theorem and the minimum image convention appropriate for periodic 



boundary conditions |Q. While the contribution from the bondlength and Lennard- 
Jones interactions to the pressure can be obtained as usual, the contribution from the 
three body bond angle interaction Vba has to be considered separately. It is found that 
the contribution of the bond angle interaction Vba can be neglected. 



V. RESULTS 



At low temperatures, T = 0.2 the sidelengths Lx,Ly fluctuate very little, and the 
system organizes itself into a configuration consisting of 12 rows and 12 columns of chain 
molecules. At T = 2.0 the sidelengths fluctuate more strongly, and these fluctuations are 
anticorrelated in time. Increases in are accompanied by decreasing Ly and vice versa 
yielding roughly constant area A. In general the measured pressure agrees well with the 
applied pressure. Small deviations are attributed to flnite size effects, and the fact that 
the pressure is measured from the virial theorem rather than from the free energy. 



Figure |T] shows snapshots of equilibrated conflgurations of 144 chains at temperatures Fig. [| 
T = 0.2,1.0,2.0 and pressures 11^ = 10,100. Each effective monomer is represented as 
a sphere of radius a/ 2. At high pressures 11^ = 100 the chains are compressed into a 
hexagonal arrangement in which all chains stand perpendicular to the substrate. At the 
low temperature T = 0.2 and high pressure the monolayer surface shows a nearly regular 
modulation resulting from extension and contraction of the chains along the molecular 
axis. At higher temperatures the surface is rough. At low pressures the chains assume a 
uniformly tilted crystalline arrangement which progressively disorders as the temperature 
is increased. At T = 2.0 the layer has molten into a fluidlike phase. 

The crystalline order and the melting of the surface layer can be seen from Figure 



^ showing the density profiles along the z direction. Figure ^ shows also results at an Fig, 
intermediate pressure Ilyi = 30. At low temperature T = 0.2 the system is crystalline. 
For = 10 (and T = 0.2) there are six pronounced maxima in the density profile. The 
position of last maximum corresponding to the tail group z = 6 is significantly shifted 
from the value Gdo = 4.2 corresponding to untilted chains. This indicates that the chains 
are tilted. At higher pressures Ua = 30 and 11^ = 100 (and low temperatures T = 0.2) the 
density profile exhibits 11 respectively 12 maxima instead of 6. This fact together with the 
positions of the maxima shows that the chains are alternatingly stretched and contracted 
along the molecular axis. The concomitant doubling of the number of monomer layers 
has also been observed in simulations at constant volume when the density becomes as 
high as 1.3 |3T| , |29[| . At the intermediate temperature T = 1.0 and 11^ = 10 the crystalline 



order is stabilized by the tilt order. For 11^ = 30 and 11^ = 100 the surface layer is 
begininning to melt. At T = 2.0 the tail group layer {i = 6) is molten for all pressures 
but the layered structure is still intact in the lower layers. Note however that for T = 2.0 
and = 10 the second layer (z = 1) is spreading into the substrate surface {z = 0) 
indicating a melting of the head group lattice. 

Defining the monolayer thickness from the largest (rightmost) inflection point along 
the density profile allows to analyze the change in monolayer thickness as a function of 
temperature and pressure. The results are collected in Table |. As expected the thickness 
increases with increasing temperature and with increasing spreading pressure. 



To analyze the tilt order we display in Figure ^ the Voronoi diagrams for the con- Fig. 
figurations shown above in Figure |l|. In these figures each chain molecule is shown as a 
dot with a short line attached to it. The dot represents the center of the head group, 
and the short line represents the projection of the head to tail vector into the substrate 
plane. In addition the Voronoi cells are drawn around each head group to visualize the 
degree of crystallinity in the head group lattice. These figures show that at high pres- 
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TABLE I. Monolayer thickness for different temperatures and pressures 





nA = 10 


Ua = 30 


Ua = 100 


T = 0.2 


3.7 


4.0 


4.1 


T = 1.0 


3.9 


4.1 


4.2 


T = 2.0 


3.9 


4.2 


4.3 



sures Ha =100 the chain molecules are untilted at all temperatures, although defects 
and randomly oriented projections of small magnitude are visible at elevated tempera- 
tures. At = 100 the head group lattice is crystalline at all temperatures. At low 
pressure and low temperature (11^ = 10, T = 0.2) the figure shows tilting in agreement 
with the density profiles. The tilt is directed towards the nearest neighbour molecule. At 
= 10,T = 1.0, however, the average tilt direction is towards next nearest neighbours. 
Both of these tilt directions have been observed in experiment [Q. At high temperatures 
(Ha = 10, T = 2.0) both the tilt order and the crystalline ordering of the head groups are 
destroyed. 

The expected magnitude of the absolute value of the tilt angle (|^|), defined as the 
angle between the end-end vector and the surface normal, is displayed over a wide pressure 



range in Figure ^ for the temperatures T = 0.2, 1.0, 2.0. At all temperatures the tilt Fig. |^ 
angle increases as the spreading pressure is lowered. However, such an increase does not 
necessarily indicate a larger collective tilt PSIE^P^PSI. An increase in tilt magnitude can 



also be caused by an increase in the fluctuations of the tilt angle. The two situations can 
be distinguished by measuring additional quantities. The projection R^y of the end-end 
vector e = {cx, Cy, e^) into the xy-plane is defined as 

R.y = VW + ^^ (7) 
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where and are the configuration average of the x and y components of the end-end 
vector e. The quantity R^y is independent of the tilt direction. Figure ^ which shows 
Rxy over the same pressure range as the tih angle, indicates that the jump in {\6\) at low 
pressures has different origins for T = 0.2 and T = 2.0. For T = 0.2 the simultaneous 
increase in Rxy indicates a tilting transition at around Ua ~ 30. The smooth behaviour 
of Rxy for T = 2.0 indicates that the jump in {\6\) at this temperature is not caused 
by tilting but by melting. This interpretation is suggested also by the analysis of the 
configurations and density profiles shown in Figures It is further corroborated by 
measuring the orientational correlations between neighbouring chains. The orientational 
correlations are measured through the quantity 



where ^atat is the angle between the end-end vectors of two nearest neighbour chains, and 
the sum extends over all nearest neighbour chains of a given chain. The nearest neighbour 
chains are defined as those six chains whose head groups have the smallest distances from 
the given chain. Because of steric hindrances the correlation K^n is usually very strong. 
Figure ^ shows the pressure dependence of K^n at T = 2.0 which exhibits a sharp drop 
at low pressures. The same curves for T = 0.2 and T = 1.0 are indistinguishable from 
the line at Knn = 1. The sharp drop in Figure ^ indicates the presence of a melting 
transition as suggested by the jump in {\9\) at the same pressures. 

The melting transition at T = 2.0 and low pressures appears also when plotting the 
specific area LxLy/N {N = 144) or the ratio of sidelengths Lx/Ly as a function of pressure 
as shown in Figure ^ The area per molecule shows a steady increase from values around 
0.72 to 0.93 and then a sudden jump to a value close to 1.2. 

At low temperatures T = 0.2 the specific area shows a sudden increase at the critical 
pressure 11^ ^ 30 for the tilting transition. This increase is associated with a similar 





12 



increase in the ratio L^/Ly reflecting an orthorhombic distortion of tlie liead group lattice 
during the tilting transition. At high pressures the ratio L^/Ly approaches the value 
2 / for the hexagonal lattice independent of temperature. 



VI. DISCUSSION 



Our present simulation results combined with previous simulations at constant area 



|, |28|j31| , |29[| indicate that the idealized coarse grained model has at least four distinct 
phases: an ordered phase without tilt, two ordered phases with tilt towards nearest and 
next nearest neighbours respectively, and a highly disordered fiuidized phase. 

At sufficiently high pressures the chains organize into a hexagonally ordered condensed 
phase without tilt. This phase persists throughout the whole investigated temperature 
range, and it can be understood as a close packing configuration. As the pressure is 
increased the individual monomers begin to move away from the minima in the bond angle 
potentials, and this leads to stretching and contraction of the chains along their molecular 
axis resulting in the visible periodic modulation of the free surface. The modulation of 
the surface is expected to depend on the relative strength of the bond length potential as 
compared to the Lennard- Jones and the bond angle potential. 

At lower pressures the chains tilt either in the direction of the nearest neighbour chain, 
or in the direction of the next nearest neighbour chain. The tilting may be understood as 
a result of the Lennard- Jones attraction between neighbouring chains. The tilt direction 
is the result of a complicated interplay between packing, energetic and entropic effects. At 
spreading pressure Ha = 10 and temperature T = 0.2 we observe nearest neighbour tilt 
combined with a periodically modulated stretching of the chains along their axis. This 
indicates that packing effects are dominant for this pressure. At higher temperatures, 
T = 1.0, the tilt direction changes into the energetically favoured next nearest neighbour 
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direction. A complete understanding of these effects will require not only more simulations 
at intermediate pressures and temperatures but in addition a systematic exploration of 
the space of interaction parameters which is beyond the scope of the present study. 

Finally at still higher temperatures the individual chain conformation and the head 
group lattice become disordered. The system enters a fluidized phase in which more and 
more configurations become accessible as the temperature is increased. 

A similar sequence of phases has been observed also in simulations at constant area 
when varying temperature and density ||27| , |28| , |3T1 , |29|1 . Contrary to those simulations the 



phases appear here in a much narrower temperature range. In particular the fluidized 
phase appears already at much lower temperatures than in simulations at constant area. 
If the rough estimate a ~ SAfrom section |T|is used to map the model to experiment then 
the area per molecule increases from values around 23A/molecule to 30 A/per molecule 
as seen in Figure |^ This agrees roughly with the experimental values for the transition 
from liquid condensed to expanded phases, and suggests to identify the fluidized phase of 
our model with the liquid expanded phase in Langmuir layers. 

Because of the wide range of pressures studied here the pressure resolution is poor, 
and hence our results do not allow to identify the order of this transition from the present 
study. Finite size scaling analyses with higher pressure and temperature resolution are 
needed to conclusively identify the order of this transition. 

Our simulation at constant pressure using an algorithm compatible with the symme- 
try breaking allows to answer the questions, posed in the introduction, concerning the 
restructuring of the head group lattice during tilting. Figure |^ demonstrates the exis- 
tence of an orthorhombic distortion of the head group lattice during the tilting transition. 
The magnitude of this effect may also be estimated from Figure |^. This distortion is sup- 
pressed in simulations at constant area because the simulation box is rigid in that case. 
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A second effect of the rigid simulation box in constant area simulations with periodic 
boundary conditions is the appearance of negative pressures at low densities and temper- 
atures The negative pressure arises from the attractive part of the Lennard- Jones 



interactions and complicates the analysis of simulations at constant area. Our constant 
pressure simulations presented here are free from such complications. 

In summary we find that our coarse grained model exhibits untilted, tilted and fluidized 
phases similar to those observed in experiment. The model reproduces the positional 
ordering, bond orientational ordering and tilt ordering similar to that observed in the 
condensed phases of the experiment It does not allow for herringbone ordering, 

however, because of the cylindrical symmetry of the chains. Our simulations at constant 
pressure show for the first time the experimentally known orthorhombic distortion of the 
head group lattice during tilting transitions in the liquid condensed phases. We have 
also shown that the choice of ensemble (constant pressure versus constant area) plays an 
important role in small scale simulations as evidenced by the large difference in transition 
temperatures for the appearance of the fluidized phase. This will become particularly 
important when trying to make quantitative predictions. 

Of course it is clear that an idealized model such as the one studied in this paper 
can only be a small step in understanding the complexity of Langmuir monolayers, and 
improvements of the model are desirable in several directions. In particular the impor- 
tant role played by the head groups and their interactions are inadequately represented in 
our model. Similarly the cylindrical symmetry of the chains is unrealistic, as mentioned 
previously. Improving the model and simulation technique step by step, and studying the 
intermediate idealized models, it will ultimately become possible to understand quantita- 
tively the full complexity of the phase structure observed in the experiment. 



15 



ACKNOWLEDGEMENT: We are grateful to Prof.Dr. K. Binder for discussions, and 
we thank the Deutsche Forschungsgemeinschaft through its GraduiertenkoUeg "Physik 
und Chemie supramolekularer Systeme" (F.M.H.) and Norges Forskningsrad (R.H.) for 
financial support. 



16 



FIG. 1. Snapshots of equilibrated configurations containing 144 molecules at various tem- 
peratures and spreading pressures corresponding to a) T = 0.2,11^ = 10; b) T = 0.2,11^1 = 100; 
c) T = 1.0, = 10; d) T = 1.0, Ha = 100; e) T = 2.0, = 10; f ) T = 2.0, = 100. 

FIG. 2. Total monomer density profiles as a function of the distance z from the substrate. 
The densities are normalized to 6. The head group monomer (i = 0) at z = is not included in 
the plot. 

FIG. 3. Voronoi diagrams of equilibrated configurations shown in Figure I at various 
temperatures spreading pressures corresponding to a) T = 0.2, 11^ = 10; b) T = 0.2, 11^ = 100; 
c) T = 1.0, = 10; d) T = 1.0, Ha = 100; e) T = 2.0, Ha = 10; f) T = 2.0, Ha = 100. The 
corresponding box dimensions are ^^^(T = 0.2, IIa = 10) = 11.87, iy(T = 0.2, IIa = 10) = 9.82, 
L^(T = 0.2, Ha = 100) = 10.92, Ly(T = 0.2, Ha = 100) = 9.45, L^{T = 1.0, Ha = 10) = 11.58, 
Ly{T = 1.0, Ha = 10) = 10.21, L^{T = 1.0, Ha = 100) = 11.15, Ly{T = 1.0, Ha = 100) = 9.63, 
L^{T = 2.0, Ha = 10) = 14.49, Ly{T = 2.0, Ha = 10) = 11.61, L^.(r = 2.0, Ha = 100) = 11.38 
and Ly{T = 2.0, Ha = 100) = 9.93. 

FIG. 4. Average tilt angle {\9\) as function of spreading pressure for temperatures 
r = 0.2, 1.0,2.0 

FIG. 5. Averaged projection of the end-end vector Rxy into the substrate plane as function 
of spreading pressure for temperatures T = 0.2, 1.0, 2.0 

FIG. 6. Orientational correlation between neighbouring chains as function of spreading 
pressure for T = 2.0 

FIG. 7. Averaged specific area L^Ly/N and sidelength ratio L^/Ly as function of spreading 
pressureforr = 0.2, 1.0,2.0 
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